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We study the renormalization group behaviour of the topological charge density 
in full QCD on the lattice. We propose a way of extracting the necessary renormal- 
ization functions from Monte Carlo simulations. 



One important issue which has come up in recent years, in the study of strong interac- 
tions, is the so called proton spin crisis The issue concerns the on-shell nucleon matrix 
element of the singlet axial current 

(P, e| jI Ip', e' ) = u{P, e) [^1(^)7^75 - ^2(^2)^^75] u{p', e') (1) 

where e, e' label the helicity states and k is the momentum transfer. In a naive wave 
function picture Gi(0) can be interpreted as the fraction of the nucleon spin carried by the 
quarks. Experimental determinations lead to an unexpectedly small value of 6*1(0), calling 
for a theoretical explanation in the context of a nonperturbative study in QCD. Gi(0) can 
be extracted from the nucleon matrix element of the topological charge density q{x) by 
using the fact that in the chiral limit the anomalous divergence of is proportional to the 
topological charge density. In the presence of dynamical fermions, one has 

(p,e|g|p',e') = MB{k^)u{p,e)t^,u{p' ,e'), (2) 

with 



NfB{0) = Gi(0) . (3) 

Preliminary lattice studies of the on-shell matrix element of q{x), in quenched QCD, have 
been done in Refs. @-§[. An unquenched study using a geometrical definition of topological 
charge density has been reported in Ref . . We also would like to note the direct calculations 
of the on-shell proton matrix element performed in quenched QCD reported in Refs. 

Notice that the above discussion overlooks the non-invariance under the renormalization 
group of the operators involved; indeed, unlike pure gauge theory, in full QCD the topological 
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charge density mixes under renormalization with d^j^. Therefore the relations of the matrix 
elements of and q{x) with the experimental observables like G'i(O) cannot be strictly 
given by Eqs. (|I]) and (|^). Instead, one must look for equations involving renormalization 
group invariant quantities. This will appear clearer in the following, when we analyze the 
renormalization of operators regularized on the lattice. 

In this letter we address the scale dependence of the matrix elements of operators entering 
the anomaly equation, in lattice regularization, and provide the prescriptions necessary to 
undertake a lattice calculation of 6*1(0) in full QCD. 

Let us first discuss this issue in the case of dimensional regularization with minimal 
subtraction. In euclidean space 

d^jlix) = t2Nfq{x)+t2p{x) , (4) 

where 

/=i 

qi^) = ^^f;,{x)f;,{x) , 

p{x) = J2^f'^f(^h5'iPf{x). (5) 

The nonrenormalizability property of the anomaly in the MS scheme means that the anomaly 
equation should take exactly the same form in terms of bare or renormalized quantities. 
However the renormalization of the three pseudoscalar gauge invariant operators d^j^Ax), 



q{x) and p{x) turns out to be nontrivial; in the MS scheme it reads |lTO|JTl 
fi2Nfq{x)\^ (I z-l 0\ /i2A^/g(x''^ 




y i2p{x 




(6) 



where, to order g^, 



167r4 8 ' e 



7 + O (,^) , (7) 



e = d/2 — 2 and Cp = ^^^^ = |. The structure of the renormalization matrix above assures 
the stability of the anomaly equation under renormalization. The associated anomalous 
dimension is given by 



7(^?) 0' 
^{g) 




^^(9) = f^^Z-' = I ^{g) I , (8) 



where 

lig) = /^^In^ = 72^7^ + 73^?' + , (9) 
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and 



72 



1 3c, 



Nf 



(10) 



We consider now the lattice regularization of QCD given by the Wilson action. Bare 
operators are defined as local analytic functions of elementary fields, i.e. Ufj_{x) and ip{x). 
At finite lattice spacing the regularized version of a continuum operator is not unique, indeed 
regularized versions of the same continuum operator can differ by 0{a) terms (at the classical 
level). We consider the following lattice regularizations of j^(x), 9^i^(x) and q{x): 



^p{x + ij)^^-i^Ul{x)ip{x) + il){x)-i^-i^U^{x)ip{x + /i) 



(11) 



5,L 



^^(^)-jJ^(^ 



/i) 



(12) 



(these definitions come out naturally from the lattice Ward identities associated to an axial 
vector fiavour-singlet transformation [jl2|), and 
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±4 

'24 X 327r2 ^ 



(13) 



(n^jy is the product of link variables around a plaquette) ||T3l , which has been largely used 



in quenched lattice QCD studies. The renormalization of in pure gauge theory has been 

a nonperturbative study for SU(3) was performed in 



studied perturbatively in Ref. 
Ref. i. 

The r-term in the Wilson action breaks explicitly chiral invariance, which should be 
recovered at a non-trivial value of the bare quark masses. Perturbative calculations and 
general arguments support the existence of a chiral limit, in which the axial anomaly should 
be reproduced In that limit, any mixings of A^j^'^ and with lower dimension 

operators, in particular ip'y^ip, must necessarily vanish, for the anomaly equation to hold. 
Indeed, mixing with ip'y^ip should be completely absorbed in a redefinition of renormalized 
quark masses; thus, in the chiral limit (i.e. 0), only the operators A^j^'^ and q^ need 

be considered in a renormalization analysis. In the following we study the renormalization 
of A^j^'^ and q^ in the chiral limit. A more general analysis for non-zero renormalized mass 
values (i.e. ^ 0) requires the introduction of a regularized version of mip'y^il:, but this 
should not lead to more than few straightforward changes in the following considerations. 

Since we are interested in the matrix elements of operators in the MS scheme, we must 
renormalize our bare lattice operators j^'^ and q^ in such a way that their renormalized 
correlation functions coincide with those obtained in MS. 
lattice renormalization matrix Z^: 



This should be achieved by a 



i2Nfq 



R 



■'tptp 



i2Nf q^ 



(14) 



It is worth mentioning at this point that, unlike the case of fiavor-nonsinglet currents, a 
nonperturbative evaluation of the matrix elements of Zi cannot be obtained by use of Ward 



identities according to the method developed in Ref. [15 . 
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Since the anomalous dimension matrix 7^(5') must depend only on the renormalization 
scheme, and not on the type of regularization, we must have 



dZL I 
^^^^ 



(15) 



This provides some constraints on the form of Z^, in particular on its logarithms. We should 
have that: (i) Zqq is a finite function of the bare lattice coupling qq; (ii) The logarithms 
appear at order in Zg^ and Z^^; (Hi) Z^q = |T^. Therefore 



^qq 



9o 



Zqi, = 90 (42 In + 4?) + <^ [do] 

Z^^ = l + gl^ll + (cJ5 In a/i + c/JJ) + O (gfj . 



(16) 



Some of the above coefficients are known. It is easy to see that S-^^ is equal to the corre- 
sponding coefficient of the multiplicative renormalization of in a pure gauge theory P^JT7[] 

^ ^ (17) 



where Zq = 0.15493. ~ 0.9084 for N = 3. From Eq. (|15D 

1 3c, 



c = c. 



(2) 



.(2) 



IGtt^ 8 



^ N 



f ■ 



(18) 



depends on the Wilson parameter r; for r = 1, we find, using the results of Ref. [|18 



/(I) 



-0.073. 



In the following we set 



R 



'i2Nfq 

Under a change of scale, the matrix elements of Oij(/x) satisfy 



(19) 



^^T^^On) = 7^ (On) 



whose formal solution is 

{0 am 



Pexp 



g(t^') 7 (q) 



(Oi?)(/x) 



(20) 



(21) 



where l3—{g) is the /5-function in the MS scheme and the symbol V indicates path ordering. 
Given the particular form of the matrix 7_ in our case, it is easy to see that V may be 
dropped. Proceeding as in Ref. [0, we choose jj,' = 1/a and write 



(0^)(l/a) = ZL{l/a){OL) 



exp 



'^(l/a) 7 iq) 



exp 



'y—ig) 



(Oij)(/i) , (22) 



and therefore 

(Ol) = 



/•9{i/«) 7 (q) 1 r fO 7 (q) 



° 7— (fl') 



(23) 



= Zf'(/i)(7o, where ^^(/i) is the couphng renormahzed in MS pO|,pT 



Notice that the vector 

M{^J){OR){^^) = exp 



° 7— (5') 



(24) 



is renormahzation group invariant, that is independent of /i, and it is what can be naturally 
extracted also from physical experiments. The lattice dependence is restricted to 



7 iq) 

L{go) ^ ZL{l/a)-' exp ^^^^ 



(25) 



The matrix M{fi) can be expanded in powers of g{fi)'. 



1 ^9if^) 



,0 ^ + W,9 

where bo is the first coefficient of the /3-function: 

1 



(26) 



bo 
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Notice that 



(27) 



(28) 



is a small number: ~ —0.0486 for = 3 and Nj = 4. 
L(fifo) can be expanded in powers of go; to order gQ : 



c „2 



1 d-^qqdo 2feo^O 

. i-(4V2f;;)^o 



(29) 



Use of the anomaly equation for the renormahzed operators simplifies considerably the 
above formulae; indeed, in the chiral limit. 



i2Nfq{ 



XI 



implying 



M(/i)(0^)(/i) 



({T)\ 



(30) 



(31) 



where 
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(T) M9w;W«)exp£^ Jgyds (32) 

We recall that 7(5') has been defined in Eq. (§). Thus, measuring the matrix element (g^) 
by Monte Carlo simulation, we obtain 

^2Njq^) = z,{gl) (T) , (33) 

where 

z,{gl) = Ln + Li2 . (34) 

Notice that (T) is renormalization group invariant. 

From Eq. (^) we find Zq{gl) = 1 - 0.8598^g + 0{g^) for = 3 and Nf = 4. We recall 
that the typical values of g^ where simulations on today's supercomputer can be performed 
are actually not small: g^ 1. Therefore as a consequence of the large 0{gQ) coefficient, 
more terms turn out to be necessary to obtain an estimate of Zq{gQ ~ 1) from perturbation 
theory. Of course an improved lattice topological charge density operator having a smaller 
0((7q) term in Zg^g^) would be welcome. The evaluation of the 0((7q) term in Zg^g^) needs the 
calculation of the coefficients d'^g^ and d^g^ in Eqs. (|16|) and 73 in Eq. (H). d^g^ requires a two 
loop lattice calculation. However an estimate may be obtained from the nonperturbative 
determination of the multiplicative renormalization of in quenched QCD, performed in 
Ref. Q, supplemented with relatively easier analytical two loop calculations of the only 
graphs containing fermion lines among those contributing to the two-point gluonic function 
with an insertion of q^. The evaluation of the other coefficients should be simpler, dg^ can be 
obtained from a one loop calculation, while 73 requires a two loop calculation in dimensional 
regularization. 

Given the difficulties of a perturbative evaluation, a nonperturbative estimate of Zglg^) 
would be welcome. To this purpose, notice that the result (|33D is independent of the par- 
ticular matrix element. We may therefore eliminate the lattice factor Zglg^), masking the 
physical signal, by measuring two different matrix elements of qi- For example, beside the 
on-shell nucleon matrix element, we may also extract the matrix element ( 0| g'^ |?7' ) from a 
lattice calculation. The wall- wall correlation function ^"(t) constructed from ( q^{x)q^{y) ) 
should have the following long distance behavior 

Furthermore 

(0|i2iV/g^|r/') = Zg{0\T\r]'), (36) 

and 

(0|T|V) = ^fml,f,>. (37) 



In the large- limit = /^r A more precise phenomenological estimate gives /, 

0.8/, mm. 



V — 
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We have thus seen how a Monte Carlo measurement of the proton matrix element of q , 
in the presence of dynamical fermions, can be related to the physical, renormalization group 
invariant quantity relevant for the spin crisis problem; the renormalization function involved 
in this process can be determined nonperturbatively from measurements of ( q^{x)q^{y) ) and 
an estimate of /^/. 
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